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ABSTRACT 


A target must choose a path between some origin and 
destination. The total travel time and the target speed are 
specified, and the target wishes to maximize the “randomness” 
Of its track subject to the spatial and temporal constraints. 
Measures of effectiveness are developed against which the 
“randomness"* of any path-producing method can reasonably be 
judged. Previous investigations into the scenario are 
reviewed and two models are developed, one using a random 
tour with drift and the other derived from Brownian motion. 
Statistics generated by Monte Carlo simulations for both 
models are compared. While the Brownian motion derivec 
process 1s not always under perfect control of the 
constraints, if the timed arrival constraint may be slightly 
violated then that process performs better against the 
measures of effectiveness and is easier for a target to 


execute than is the random tour with drift. 
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I. INTRODUCTION 


This thesis will investigate two methods for constructing 
random target tracks between two specified endpoints. The 
target is constrained to begin its journey at one of the 
points and end it at the other after a stated elapsed time. 
Furthermore, either the mean target speed or the maximum and 
minimum speeds will be specified. 


Con¢eicer the extreme case in which the distance between 


the two endpoints, the target maximum = speed, and the 
specitied elapsed time are such that the target is 
constrained to travel directly to its destination. Here 


there can oe no randomness to the path, and an attacker needs 


very little target position information in order to make 


perfect oredictions about future target position. Now 
imagine a less constrained situation in which there 15 
considerably excess time. The target is able to travel 
arcuna randomly during most of the period. It might choose 


to move quite randomly at the beginning, suddenly realizing 
at some critical point that it has just enough time remaining 
to travel directly to its destination. However, it 1S more 
intuitively sensible in this controlled time of arrival 
scenario to "spread" the randomness and effects of the 


constraints evenly across the time period. 





The two scenarias above give insight about what 
characteristics of target travel are important. In the 
completely constrained case, the target’s position does not 
vary at all from its expected position, resulting in = an 
easily inferred position throughout the scenario. 
Furthermore, the target always points direcly toward its 
destination, as well as directly away from its Seen For a 
target that wishes to Keep its origin and destination secret 
for some tactical! reason, compietely constrained travel is 3 
giveaway. On the other hand, 2 target that travels about in 
an entirely random manner may rever get to its destination. 

It 1s clear that constrained target travel between two 
endpoints invoives many tradeoffs among various constraints 
and choices. This thesis will “develop ai notion of 
"randomness" by devising measures of effectiveness that 
logically follow from the scenario and that also have strong 
intuitive appeal. Some previous apsroaches to the problem 
will be discussed briefly and then two new “recipes” “for 
target travel will be developed in detail and evaluated 


against the measures of etfectiveness adopted. 
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TI. MEASURES OF EFFECTIVENESS 





The previous scenario describing totally constrained 


target travel made apparent two dangers to the evasive 


target. First, the target risks divulging both its origin 
and destination because it always points toward ts 
destination and away from its origin. It follows that the 


less constrained a target is, the less it necessarily coints 
to destination and away from origin. Exactly how much course 
freedom is gained as the target becomes less constrained will 
be investigated in detail later. One might ccunter the 
importance of this pointing by claiming that, at best, the 
attacKer obtains only a line of bearing to the origin or 
destination. While true, this argument neglects the 
possibility that other targets may pass the same way, with 
the same origin or destination. Sooner or later the attacker 
Will get lines of bearing that cross with a ;requiarity 
sufficient to specify the critical positions. Since origins 
and destinations might be important enough to be Kent secret, 
one reasonable measure of effectiveness against which to 
evaluate any set of target paths 1s the absolute angle 
between the target’s present course and the course from the 
target’s position to its destination sampled at specified 
time intervals (c.f., Figure tia). For the totally 


constrained target this angle will always be zero. In order 
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to compare two different path generation procedures it will 
be helpful to preduce a plot of this measure against time if, 
in fact, it is time dependent. Otherwise, it will be useful 
to calculate the mean pointing angle, its standard deviation, 
and perhaps an empirical density function. 

The second undesirable quality of a completely 
constrained path is that the target’s position always falls 
exactly on its baseline position, where baseline position 
(PY is defined as the point where the target would be at any 
time of its Journey if it travelled directly to its 
destinatian at a constant speed. For example, if the distance 
between origin and destination is sixty distance units and 
the <¢pecified elapsed time for the problem is thirty time 
units, then the rate of travel along the baseline is two 
distance units per time unit. Ten time units after the start 
of the journey, the target’s position is on the baseline 
twenty distance units from the origin and forty distance 
units from the destination. The baseline 1S important 
because the more constrained the target 18, the closer it 
must <etay on the average to the baseline. Increasing 
constraint on the target by decreasing allotted travel time 
or by increasing baseline length straightens out the target 
path with the result that inferred positions are easier to 
obtain. Very importantly, even if the target 1s not forced 
to stay near the baseline as a result of being constrained, 


but chooses to do so by hovering or zig-zagging along the 
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line, an inferred course and position are easily obtained. 
Even worse, such an interred course will most likely point to 
the destination or from the origin. Acrordingly, the second 
measure of effectiveness shall be the mean square radial 
distance Between the present tarcet position and the baseline 


target position at specified time intervals during a journey 


2 
¢ J]. 


The squared radial distance has been chosen rather than the 


(c.f., Figure 1b). This statistic will be denoted as EIR 


radial distance because previous investigations into the 
problem, which will be cited in the next chapter, have tended 
to use radial distance squarad. It 1s easier to handie 


analytically than i¢ R and can be viewed as the sum of the 


t? 
x and y component squared distantres from the component 


baseline positions at any time t: -« 


It would be not difficu!t to integrate R,* over the whole 
time interval for an individuai path in order to get a scalar 
value for this performance measure. Likewise, the mean 
square radial distance ¢ ECR, “1 ») could be integrated for 
each particular path-producing recipe in order to obtain a 
scalar measure for that method. Then, two different path 
methods might be compared and, all other measures being 


equal, the method with the highest value selected. However, 


the reduction of this measure of effectiveness to one scalar 
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can result in great loss of information. Two path generation 
procedures might have the same value by the scalar measure, 
but Be characteristically different. One might tend to fave 
relatively large values of EER, “1 near the origin and 


destination, while the other method might have small values 


there, with large values for t in the middle of the travel 


period. Furthermore, each of these two different sath 
recipes might be more desirable than the other under 
different circumstances. For instance, hit) ¢ 5). 5 ceasst 


important to Keep secret the origin and destination becauéze 
they are already Known to the enemy, but it is important to 
be as undetectable as possible (or if detected, as difficult 
as possible to redect>) in between origin and destination, 
then the target would probably prefer the latter cf the tuo 
path-producing methods. In any case, for any path-preduciang 
recipe, the captain of the target should be able to !s0k at 
the plot of ELR,“} as tunction of time. Thus, the measure of 
effectiveness itself shall be either a plot of ELR,~3 against 
time, the function that describes it as such, or 4 iist of 
ordered pairs with time as the first element and mean square 


radial distance as the second. 


A careful observer might object to this measure 093of 


effectiveness. While it is true that a highly constrained 

target will exhibit a smal! ELR, “2 for all t and have easily 

inferred positions, this does not necessarily mean that 2 
Z 


large ECR, J] will guarantee that target position will be 


iS 





difficult to infer. In fact, if maximizing ELR, “1 over the 
entire time period is taken as the only measure of 
effectiveness, it should be almost as easy to infer future 
position from a few past positions as it is in the totally 
constrained case. Imagine a baseline ten distance units 
long, a maximum target speed of four distance units per time 
unit (4 d/t) and a snecified travel time of four time units. 
It seems sensible that the target would travel at maximum 
speed the entire time in order stay as far away from its 
expected position as possible. For a good analogy, imagine a 
stick ten distance units long with a piece of string sixteen 
distance units long attached to it, each end to each end. 
Then, the target might choose a track to maximze the area 
between the stick and the string. (The analogy is not 
perfect however; in the sticK and the string example the 
radial distance (Ry? integrated over the entire time period 
is being maximized rather than the radial distance squared 
cR,* 
straight away from the baseline at about a fifty degree angle 


>. One possible path choice is to drive the target 


for eight distance units, then turn back and go directly to 
destination. The area of the isosceles triangle thus formed 
with the stick as the base is approximately 31 distance units 
squared. But a clever hard working target could drive in a 
circular arc, thus enclosing approximately 39 distance units 
squared. This path, while producing a large ELR, “1 at each 


time ¢, 1s not very “random” and might be easily targeted. 
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The third and tinal measure of effectiveness is important 
in at least two ways. The first pertains directly to the 
shortcoming of the second measure of effectiveness taken by 
itself. Given a certain recipe for producing paths, it is 
necessary to produce many of them in order to obtain 
statistical estimations of both of the two measures of 
effectiveness already chosen. It may be possible to apply 
some non-linear regression method to any path, and to extract 
an estimate of future position to use as the expected 
position for any time t, rather than use baseline position as 
already defined. Such a method would quickly identify the 
smooth curves produced by any path recipe which maximized the 
mean square radial distance only. But the task would be quite 
difficult and costiy, especially when applied to the several 
hundred oaths necessary to produce good statistics. In order 
to save effort anc money, the third measure of effectiveness 
will subjectively judge how representative paths from each 
generating method "look". Against this measure, any path 
which maximizes EER, “3 alone will be rejected at a glance if 
it exhibits long straight legs or a predictably curved path. 
However, if a path does not look so regular as to be 
predictable, then a large EIR, “1 is desirable because it 
means the target is staying away from the critical line 
between origin and destination. 


A second very important reason that paths will be 


visually inspected is that they must be able to be executed 


Ze 





by an actual target. Hence, an experienced person must judge 
whether or not a path is practical regardless of its 
statistics. For example, a high-curvature path which 
maximizes ECR, “1 may not be implementable. 

There are, then, three measures of effectiveness that 
shall be used to judge path-producing recipes. Though they 
were introduced in a different order to facilitate logical 
development of them, they will be applied to a path procedure 
as follows. First, a representative path must “look good" by 
being practical, executable, and random looking (no 
noticeable regularities). If a path can pass this first 
important test, then several hundred will be generated using 
the same procedure, and the mean angle ‘and its standard 
deviation) between present couréee ant course from present 
position to destination will be calcujiated at specified time 
intervals. Then, the mean square radial distance between 
present position and baseline positicn will be calculated at 
each specified time interval. The ECR,“ wil] thus be 
estimated at various stages of path cempletion. I¢ two path 
generating recipes both produce paths that meet the origin, 
destination, time, and speed constraints, and "look good" by 
the first measure of effectiveness, then the path which 
exhibits the least pointing to destination as judged by the 
second measure of effectiveness and the greatest EIR,“ at a 


given stage of path completion shall be judged to be the more 


desirable path. Some comparisons will no doubt result in a 
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Situation in which both procedures produce paths which *1lo0o0k 
good", but each of the two methods has a better evaiuatian 
than the other in one of the other two measures. One might 
try to form some weighted combination of the two measures, 
but this is dangerous because measures of effectiveness do 
not combine well; at minimum, their units are not generally 
on the same interval scales. Such reductionism is aot 
necessary anyway. One need only regard both of the last tua 
measures and decide when it is advantageous to weight one 
subjectively over the other. In fact, any formal weighting 
system would probably not be able to capture al! scenaric 


dependencies as well as a subjective weighting. 
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III. PREVIOUS INVESTIGATIONS 


While several approaches to the problem of randomizing 
target motion are possible, one of the most appealing is 
investigated in detail by Washburn (Ref. 1]. In his model a 
target takes a random tour by choosing its direction of 
travel from a unitorm (8,20 probability distribution and its 
length of travel on the selected heading from an exponential 
distribution with parameter X (mean number of turns per time 
unit). Thus, turning points are the jumps of a Poisson 
process with parameter i, and at any point during the process 
both the backward and forward recurrence times are themselves 
distributed exponentially with parameter xX. This property of 
"“memorylessness"” is very appealing. The probability that the 
target does not turn by time (ttad), given that it has not 
turned by time t is the same as the probability that the 
target will not turn during the time interval (8,a). In 
other words, in deciding at any point in time when to turn in 
the future, the target does not remember how long it has 
already traveled on the same course. Hence, an attacker may 
not infer either that the target will stay on a new course if 
it just turned to it, or that the target will turn soon 
because it has been on the same course for a long time. 
Because of the desirability of giving an attacKer so little 


information, each leg of target travel is chosen 5y draws 
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from the following distributions: 


1/2 69 £<¢9 $ 24 
#«(O> = 
6 otherwise; 
C2) 
Neon t > 8 
€€%) = 
6 otherwise. 


After much rathematical manipulation, Washburn derives 
the probabilty dens: ty of the target’s radial distance from 


its origin, given no initial course information, to be: 


2 


$Cr, CHL 1/2) “I LAt/¢ 1-7 * *expl-avatci-(t-r4) 1. (3) 


He also makes soaine iriteresting observations: 


1. The larger xX is the more the distribution piles up 
around the origin. 

2. Given no turns, the target is uniformly distributed on 
a circie of radius Vt (throughout this thesis, a ‘V’ 
will denote a scalar which is the magnitude of the 

~~ 
vector denoted by ’V%). 

3. Very remarkably, Qiven that two steps have been 
completed at time t, the target is uniformly 
distributed in the circle of radius Vt, not including 


the circumference. 


The first two observations go hand in hand. Recall from 


Chapter I that a target which hovers about the baseiine or 
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Zig-zags along it has a smal} ECR, 73 Since it is rarely far 
+rom its baseline position. Washburn’s distribution 
quantifies exactly how much the target remains near its 
expected position relative to how much the target "hovers", 
as indicated by the parameter lambda. The only difference is 
that Washburn’s model has no drift, resulting in the expected 
position being a single point rather than a baseline of many 
points as in the drift case. A turning parameter of lambda 
equal to infinity creates the degenerate case in which the 
target is distributed on top of its oricin with probability 
equal to one. The other degenerate case occurs in 


observation two when the target does not turn and 1s then 


necessarily on the circie of radius Vt. Just as noted in 
Chapter I, BORE 3 may be high im this form o+ degeneracy, but 


given any two bits of position information, future position 
can be inferred perfectly because the ‘target has been 
travelling in a straight line. 

Though it is interesting that the target is uniformly 
distributed on the disk of radius Vt if it just happens to 
have finished leg two by time t, this fact really does not 
help a target evade an attacker. Given the set of 
instructions by which paths are generated, it would be purely 
coincidental that the second leg is completed right when’ the 
attacker looked for the target (the probability is zero). 
One might say that the target should plan to finish the 


second leg just when the attacker is expected. For example, 


me 





the target might Know when it conducts some type of svoluticn 
that makes itself more detectable. Also, it aiso mignt Know 
from how far away the attacKer must come and at what speed, 
so it has a good idea when its window of vulnerability is. 
But if it uses this information to plan two legs then it is 
not following a random tour as prescribed by the probability 
distributions set forth. If the target is going to break the 
rules, it may as well just pick a point at random inside the 
circle of radius Vt and head for it. This procedure wil! 
QUuarantee that the target has sampled uniformly over the 
disk. 

BelKin (CRef. 2] does pioneering work in comparing 
Washburn’s random tour process with a Gauss-Markov process 
for describing diffusion, the Ornstein-Uhl enbeck Ci OU) 
process. As a result, BelKin finds that the mean square 
radial distance of a target from its expected position at 


time t in a random tour process is: 


EIR, “1 = 200°/*) Can text Le. C4) 


Belkin [CRef. 3] further embellishes his analytic work on the 
random tour process by deriving the mean square radial 
distance for a random tour with arbitrary course chance 


distribution to be: 
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Here FC(@) is the cumulative distribution function of 9 for 
which there may or may or not be a proper density function. 
The subscript w on wae indicates that the velocity is total 
Velocity through the water, comprised of both a drift 
component and a randomizing component which will be explained 
later. . Notice that for the uniform distribution of the 
Washburn’s random tour model, ae = Ny = @, and equation ¢3) 
results as expected. 

The notion of an arbitrary course change distribution is 
important because the expected position of the pure random 
tour process is the origin at time zero. As such, the 
process will never cause a target to migrate toward its 
destination. Therefore, it 18s necessary to introduce some 
bias into the selection of the courses in order to Keep the 
target moving in the correct direction toward its 
destination; and even a course distribution alone will not 
cause a target to visit the destination without the addition 


of some further constraining process. While an infinite 
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number of arbitrary course selection distributions are 


possibie, severa: look promising when applied to the present 


problem. Loane (Ref. 4] suggests: 


#00) = 1/2% (1 + acose@) for @ ¢ (86,29 (7) 
© -<-04 2 1)% 

where Bias in the positive x direction can be controlled by 
the selection cf the parameter «a. Because constraining 
target travei by course distribution alone only causes the 
target to migrate toward its destination, without giving any 
guarantee that the endpoint will be visited, some process 
must be devised to meet the visitation constraint. One of 
many possible course distributions, which also injects the 


time constraint into course selection, is: 
6 —~s normal (CUS, ArcostT /T >] (8) 


where CUS is the course to the destination at the beginning 
of the leq, VT. is the time it would take to go directly to 
the destination at mean speed or maximum speed (whichever is 
selected as important for a particular rund, and - is the 
time remaining in the problem. Notice that in the most 
constrained case in which ue = Ts the variance of the course 
distribution equals zero and the target travels straight to 
its destination. While this course distribution causes the 


target to travel over the destination, it does not control 


when it will cross it. Hence, it does not solve the 


Figs. 





controlled time of arriva! scenario devised in Chapter I. 
Belkin CRef. SJ] meets the visitation constraint with the IQU 
process and is able to distribute the contraints evenly 
across the whole path through optimal control of velocity. 
This means that the target arrives at destination when it has 
to, without either traveling to it straightaway and hovering 


for the remainder of the period, or by randomizing travel 


until the last possible opportunity to get straight to the 
destination on time. Nonetheless, mean square radial 
distance from baseline position will necessarily need to 


begin decreasing at some point in order to guarantee that the 
target visits its destination. Two observations made by 
BelKin motivate the approaches taken in this thesis for 
devising path recipes. He States iKef. 6] that "...it is 
possibdle to approximate a random tour orocass with arbitrary 
course change distribution by an I10U process with Ilinear 
drift.* He also notes (Ref. 7) that "...as the constrained 
process approaches the terminal constraint at T, the process 
behaves precisely like an unconstrained process running 
bacKwards in time from T." T, in this instance, is the total 
time alloted for the problem and the terminal constraint 15 
the requirement for the target to be at the destination at T. 

As aresult of BelKin’s insights, the first of two path 
producing recipes that will be developed in this thesis will 
treat constrained travel between two endpoints as a random 


tour process with drift (distance between endpoints divided 
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by total time), executed from both endpoints with an attempt 
to connect the two separate paths somewhere in the middle. 
The drift vector applied to the path beginning at the origin 
will “Blow" the target toward the destination, while the 
drift for the path starting at the destination wil! “blow" in 
the opposite direction, toward the origin. This approach 
constitutes a discrete approximation of the IOQU process with 
drift, wherein the random tour out of the destination is one 
half of the total problem run in reverse. 

BelKkin’s solution to the problem of constrained tarqet 
travel is the IOQU process and he has developed a computer 
simulation named IQUTRK for which he presents some sampie 
paths (Ref. 8]. While the genesis leading to his adoption af 
the [0U process begins with Washburn’s random tour, Belkin 
Proceeds to attempt an approximation of the random tour 
process using Brownian diffusion (Refs. 9,10]. However, fis 
investigations lead him (URef. 11] to the very ciever 
realization that the functional form of the mean = square 
radial distance, ECR,“1, is exactly the same for the [QU 
process as for the random tour after making only two simple 
parameter substitutions (Ref. 12]. He concludes his analysis 
[Ref. 13] by stating that "...if one is constrained for 
theoretical or computational reasons to approximate the 
motion of a randomly touring target by Gaussian diffusion, 
then the Ornstein-Uhlenbeck displacement model is to be 


preferred to the Brownian motion model." 
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The random tour process is important and has many 
desirable properties that will be discussed further when the 
first new path recipe in this thesis is developed and 
analyzed. As indicated previously, the random tour will not 
b2 approximated by any other process; it will itself be 
executed in a new way. The second model to be offered in 
this thesis will examine Brownian motion, not as an 
approximation of a random tour, but as a basis for solving 
the constrained target motion problem in another way. Though 
rejected by BelKin as an approximation of the random tour 
process, Brownian motion stands alone as a method for solving 


the problem at hand in a very novel way. 
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IV. PROCEDURE 


A. RANDOM TOUR WITH DRIFT 

In order tc evaluate the random tour with drift against 
the measures of effectiveness delineated in Chapter II, a 
Monte Carlo simuiation is devised. Specifically, two random 
tour processes are axecuted, one from each endpoint, in order 
to Qguarantee that the destination Cactually the starting 


point of the secand process) is visited as required by the 


constraints. Recali from Chapter III that the reason for 
executing two sedarate processes stems from BelKin’s 
observation that th2 path approaching the destination looks 


1iKke an unconstrained random tour run in reverse. Hence, the 
path which visits the destination is constructed by executing 
an unconstrained random tour originating at the destination. 
For all paths, and without loss of generality, the origin is 
zero on the x axis of a Cartesian coordinate system in two- 
Space and the destination 1s however many distance units 
desired in the sesitive x direction. The random tour 
beginning at the origin is “bBlown”* in the positive x 
direction by the drift vector and the process originating at 
the destination is “blown” in the negative x direction by an 
opposing drift vector of the same magnitude. The goal is to 
Generate two paths that "blow" into each other somewhere 


between origin and destination. At the outset, the distance 
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between origin and destination, the time period, the mean 
target speed through the water Vs and the target turning 
rate (lambda) must all be specified. Orift velocity is then 
in the positive or negative x direction, depending on which 
of the two random tours is being executed, with a magnitude 
equal to total baseline distance divided by total problem 
time. Figure 2 below illustrates the other important 
velocity, the randomizing velocity oS the magnitude of which 


must be computed before the random tour can begin. 


Figure 2: 


Water Velocity as a Compositicn of Drift 
and Randomizing Velocities 





Because the target must have a mean speed through the water 


of Voy? and V is determined by baseline length and total! 


drift 


run time, there is no freedom in choosing mean hnue Recal] 
that the angle theta is the direction of target travel chosen 
from the uniform probabilty distribution (8,20 and notice in 
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Hence, given V4 and a desired mean an the randomizing 
velocity for the random tour with drift should be: 
US = Uys (18) 


rc Ww d 
Notice that as the drift speed approaches the water <zpeed, 
the randomizing speed goes to zero. This is reasonabie 
because when the target has just enough time at a given speed 
to go straight to the destination, all its velocity aqnes 
toward matching the drift velocity and there is no excess 
velocity left over for any randomization; the target is 
totally constrained. 
With these preliminaries taken care of, the simulation is 
executed beginning with the “left” path Cout of the origin at 


zero on the x axis) according to the following procedure: 


1. Choose the direction of travel from UC8@,2D. 
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Zz. Select the length af each leg of the process as Vit, 
where t is determined by a draw from an exponential 
distribution with parmeter lambda. 

3. Compute the actual vector v describing travel for the 
leg by adding a prorated drift vector U, the magni tude 
of which equals (Ct / total problem time) x baseline 
length], to v which was determined in steps 1 and 2. 

94. Continue executing steps 1-3 until the cumulative 
travel time for all left side legs is at least .475 of 


the total time for the whole problem. 


After one left side path is simulated, a similar process is 
executed from the "right" side (Cdestination>) using an 
opposite drift direction. As a result there are now two 
random tour paths which migrate toward each other at the 
determined drift rate. The two processes use up at least .95 
of the total time alloted for the target travel from left to 
right. Now, a circle with radius .85 x total time x desired 
water speed is drawn around the left endpoint. The 
proportion of time left in which to join the two sides is not 
arbtrary, but depends on the length of the total time ot the 
problem. Here, the total time is assumed to be between 394 
and 6@ time units, maKing the joining period between 3 and 6 
time units. During this connecting period the target should 
maneuver as suggested by X, laying down a path that will make 


the connection. If the right path endpoint falls within the 
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circle, then a “maitch" has occurred and the twa connected 
paths constitute a complete sample path satisfying al! 
constraints. If a match is not realized in a specified 
number of iteraticens of the right path, then a new left path 
is executed and another attempt to generate a matching right 
path is made. The number of right paths that are generated in 
attempting to match a given left path and the number ef left 
paths that are generated after a right path match iteration 
Jimit is reached can coth be controlled by parameters in the 
simulation. The entire process is continued until 188 
matches ‘or any other number desired) are obtained or until 
the Jimits for the number of attempts are reached. The 
motivation for this precedure of matching is twofold. First, 
desired statistics can be generated to test the symmetry of 
both sides since many right side and left side paths will be 
generated during the quest for 1866 matches. Secondly, in 


attempting to get matches, the number of right sides al) owed 


in order to match a given left side can be strictly 
controlled. The lower the limit on right sides, the more 
likely it will be necessary to generate another left side. 


In this manner, the total number of left sides generated in 
order to obtain 1846 matches can be much greater than 1866, and 
the number of right sides will be even greater yet since 
right sides are constrained to match left sides and not vice 
versa. This path generating method results in two complete 


distributions of paths for the left side. Some left hand 


33 








paths result in matching right hand paths, white other left 
paths are rejected because no right paths can be generated 
Cwithin Jimits) to match them. The retained left paths 
comprise a set which will be called the "clipped" set, while 
all paths, rejected or not, will be called the “unclipped” 
set. These names are reasonable because the paths that match 
are more constrained; they have not seen rejected as 


unacceptable. A comparison of the statistics of interest for 


both left hand sets will indicate how much degradation 
results from constraining the clipped set. Additionally, 
there is a third set of paths which contains all the right 
paths which match left paths. Statistics are generated for 


this set in order to check symmetry with the left hand 
clipped set. ° 

In order to facilitate the generation of valid statistics 
it 1S necessary to sample the left and right paths at the 
same specific time intervals. While 1% i353 much easier to 
sample a path at the end of each leg, the resulting 
statistics do not provide consistent comparisons along the 
time line. Accordingly, statistics are qenerated at times 
determined as a function of the turning rate, A, so that the 
probability of generating statistics twice on the same leg 


is .9@61. The statistics generated at each time check for al] 


three sets of paths are: 


1. Mean magni tude of the angle between present course and 
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course to destination from present position; variance 
and standard deviation of that mean angle. 

2. Mean square radial distance between baseline aosition 
(pure drift prorated for cumulative run time) and 


actual position. 


Using the procedure described here, statistics are generated 
for a water speed of 4.8 dist/time and drift rates of 3.75, 
3.8, 2.8, 1.8 dist/time, all for turn rates of 1.80, 2.8, and 


4.6 per time unit. There are a total of twelve cases. 


B. BROWNIAN-DERIVED MOTION 

Since Robert Brown first observed the highly irregular 
motion of suspended pollen particles in 1827, mathematicians 
and scientists have spent a great amount of time and effort 
investigating Brownian motion, and its widespread 
applicability to naturally occurring events makes if a 
reasonable candidate for randomizing target motion. Prom a 
probabalistic viewpoint Freedman [Ref. 14] defines normalized 
Brownian motion to be a stochastic process (BCt):8 §¢ t < %} 


on a sample space N with properties (a-c) for points 8 ¢ : 


Ca> BC8,W) = 8 for each QO, 
(b> BC.,m is continuous for each WW, 


(cd) for @< t, < t., Sone be ass c ti the increments B<t >; 


Bae) — B(t >, -->, Bct) = BCt 


> > are independent and 


t 
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normally distributed, with means @ and variances tis 
t, 7 Cpoeees La a Co oar 


Additionally, this process is Markovian, meaning that the 
future is conditionally independent of the past, given the 


present; 
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where the 1’s are elements of the state space for the random 
variable xX. 

(ince again, the problem of constrained target motion 
shall be framed in two-space using a standard Cartesian 
coorclinate system. Taking position to be specified by each 
component separately, the conditional probabilities for each 


coordinate at time t, given the initial constraints are: 


P£ X¢t) ' X(@) = 8, X(T) = L? 


il 
x 


Giz 


Pl YCt) =ys YCO @, YCT) > 


where t is cumulative run time, T 1s total alloted time and 
wes t $$ TD, L is the positive x coordinate of the 
destination, and 6 is the y coordinate of destination. 
Hence, the Cartesian layout is exactly as it is tor the 
random tour case. Conditioning the Brownian motion 
constructs a “Brownian bridge” between the origin and the 
destination, by which the initial problem constraints are 


met. The x component conditional probability expands to: 
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fase =; ACD = 6 } ser C KCTS KC t) = x 3 
P€ X(T) =bL i X(@) = @ } 
SES9, 
=P i X(t) Ss x, eXtT = L i xc = 9 > 
P € X(T) = Li X(@) = @ 3 


and the conditional density when X(t) is distributed Gaussian 


becomes: 
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where: 
2 
Q= (x*/0*) + EtL=x> “7We" (Tt) 13 = (L276 I>. 


After completing the square for @ and simplifying, the 


conditional density reduces to: 
ot ; a 
exp€—1/20 Cx-Lt/to /¢ oe tC 1-t/T9) J 3 


cop" 2 oo . etli-ctyT) 313". 
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Hence, X¢t) and YCt) (by a similar derivation] are distributed: 


XC t} -~ Normal C€ Lt/T, otc = t/Toe 


Ve NormalmionmoSticie=seaT) 1. 


€ 16) 
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In both distributions of is a Dhysic2l parameter that may be 
specified, and the total variance describes a pezrabola when 
plotted against t. The maximum variance is 7). (7.74) and 
occurs at t = T/2. Very importantly, the variance is zero at 
both the origin ¢ t = @) and the destination ( + =T>). 

One of the measures of effectiveness chosen in this 
thesis against which to judge a Sn path exanines how a 
path "looks" and whether it is executatie by ai crew. 
Brownian motion, as defined, is not executable because it is 
a continuous process with an everchanging velerity, and as 
such, 1s impossible for a large target to duplicate. However, 
there 1S no reason why the continuous process cannot be 
sampled at various times and the selected points be made the 
endpoints for leqs of sfraight line travel. Of course, as the 
time interval between samples 1s lengthened, the linear 
approximation connecting the sample pcints hecomes less 
"Brownian", Nonetheless, the process sti!!] retains vestiges 
of its Gaussian properties, and the fiavor of Srownian 
motion. 

An exponential distribution should te used to determine 
the sample times for the same reason that is advanced for the 
random tour: the memoryless property applied to course change 
1S advantageous to the target. Exactly how to sample from 
the normal distributions (16) 1s the last problem to be 
solved before a procedure can be devised for executing paths. 


One might sample by generating a string of exponential 
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random variables using the desired parameter lambda in order 
to get the time jumps ¢ tis Cases Ce > of the Poisson 
process which determines course change; then make a draw from 
each from the x and y normal distributions at times t = ti 
ty + toseees ae + al to obtain x and y coordinates for the 
path legs. While the mean of the y coordinate draw is always 
zero, there is nothing in this method to prevent cne_ y 
coordinate draw from falling on one side of the x axis and 
the mext draw on the other side. As the cumulative time 
approachs T/2 the variance of the draw becomes relatively 
large and it is quite likely that two consecutive draws which 
fall on opposite sides of the mean will be very far apart 
relative to the corresponding time step. A similar arqument 
applies to the x coordinate, and because position is. 
determined by Both coordinates, if two consecutive draws are 
on the opposite sides of their means for both the x and -y 
coordinates, the distance between two consecutive positions 
could be very far apart. The result is that a ridiculsusiy 
high velocity is required in order for the target to travel 
from one sample point to its successor in the given time 
step. What is needed is some way to guarantee that a 
successive position is tied to the one before it. Future 
position must be conditioned on present position, and that is 
what is missing in this method. Qne possible remedy might be 
to reflect each draw across its mean if needed. Thus, the 


Variance of the draw is preserved while the distance between 
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zuccessive positions is shortened to a manageable length. 
However, this type of optional reflection does not capture 
the true Markovian property of each draw, but it does provide 
insight into how each successive draw should be done. 

Suppose that the first draw is carried out without 
reftiection as described above. The first leg will connect the 
origin with a point selected by a draw from the two normal 
distributions ¢16) substituting t = tT) (the first time 
selected from the exponential distribution). InvoKing the 
Marxovian property, all information prior to this latest 
point is superfluous. Hence, there exists a completely new 
controlled time of arrival problem, wherein the latest point 
becomes the new origin and the line between that point and 
the destination becomes the new baseline. The total run time 


for the new problem is T -t This) iterative process 


ir 
continues until the original total time ¢€T) has expired. On 
the final leg, the cumulative time for the two coordinate 
draws 18s t = t, 5 t. feet oa = T, so that the variance in 


C14) goes to zero and the target visits its destination = on 
time. This procedure 1s a sensible one that seems to capture 
the properties of Brownian motion while at the same time 
producing discrete linear legs of target travel. Thus, it is 
the procedure adopted to produce the second type of path 


developed in this thesis, and the steps of the recipe are: 


{. Select baseline length (distance between origin = and 
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destination = LL), total run time (T), turning rate 6d, 
maximum and minimum target speed, and Physical 
Parameter for variance co Simulation frame of 
reference is a Cartesian coordinate system in 2-space, 
with origin at (8,6) and destination at (L,@). 

Generate a string of random variables from an 
exponential distribution with parameter lambda, and 
truncate the nth vaiue in the string so that the sum of 
the random variables t + t., +eoot te = T. 

Make a draw from the normal distributions in (16), 
using t = cumulative run time. For the first draw, 
t = t 


for the second draw, t = t, + and so on. 


weer: 
Measure distance between new point (Xi s1;? and previcus 


point CK >. If the distance/t, 15 greater than 


Ve 
the maximum target speed, truncate the leg as 
illustrated in Figure 3 if the distance/t. is less than 
the minimum target speed, extend the leq. Either 
truncation or extention results in new (X57). 

Reframe problem as a completely new one, using CX. 515? 
as the oricin, (L,@% as the destination, and T= time 
remaining. Execute steps 3,4 again. The last draw 
forces (X,Y) to be equal to (L,8) unless truncation or 
extention occurs, in which case problem ends on time 
with target short of destination, or target continues 


to destination, in which case total travel time exceeds 


T. One path is completed. 
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Figure 3: 


Leg Truncation and Extension 


¥ Figure 3a: Truncation 
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6. Generate several hundred such paths and compute the 
statistics of interest for each one at specific time 
intervals ‘(mean square radial distance between present 
position and baseline position of original problem; and 
absolute difference between target present course and 
course to destination from present position). Look at 
plots of paths to determine if they "look good", as 


discussed in Chapter II. 


After statistics and plots are generated for both the 
random tour and Brownian motion derived paths, the task is to 
compare the two approaches in order to determine their 


respective strengths and weaknesses. 
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V. RESULTS 


&. RANDOM TOUR WITH DRIFT 

One measure of effectiveness for judging paths requires 
that they “look good", as detailed in Chapter II. It makes 
sense to try this measure on a path procedure first before 
expending effort to evaluate amass of statistics. If the 
candidate paths can be rejected on sight then time will not 
be wastec on the other two measures. 

Figures 4a-f are representative paths generated by the 
random tour with drift method. The circle in these figures, 
as described earlier, is inscribed with the end of the left 
path as its center and has a radius = Nae x time remaining in 
probiem after the execution of both the left and right paths. 
Thus, the target can easily travel from the end of the left 


faq to the end of the right leg in the allotted time at the 


a 


tated water speed. The ratio MA listed on each of the 
ti:gures measures just how constrained each path 15s. The 
ratios ranges from zero for the unconstrained case (random 
zcur without drift) to positive infinity for the totally 
constrained case (straight line between origin and 
destination). Figure 4a shows a path for which eer = 2.673 
clearly this path does not "look good" and is not acceptable. 
It is almost a straight line between origin and destination, 


and future target position is easily inferred from only a few 
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fixes. Of course, any method for devising a path under such 
constraints is bound to fail. A target cannot randomize its 
motion if it must travel directly to its destination. The 
rejection of the path in Figure da is not a rejection of the 
idea of random tour with drift, but only of the constraints. 
Even though 2.69 is in the low end of the range of Waewee any 
ratio greater than about 1.5@ produces paths that are 
unacceptable. In Figure 4b, where the ratio is 1.13, the 
path is still highly constrained but looks much better. The 
Grift constraint really begins to loosen in Figure 4c, and 
in Fiqures 4d-f the ratio is 6.26. Here, the randomizing 
qualities of the random tour are evident. The paths al] 
"leo0K good", and whether they are executable by a target 
Gepend:s only upon what the time and distance units are. For 
instance, if the time unit is one half hour, Figure 4f has 
the target turning every 7.5 minutes on the average. This 
turn rate is not realistic for a large target. However, the 
path in Figure 4d has a mean turn rate of two per hour and is 
certainly executable. The main concern at this point has 
been put to rest; the paths in Figures 4c-f for which the 
time constraints are reasonable, 100K good enough to warrant 
further evaluation. 

The next measure of effectiveness against which to 
evaluate the random tour with drift is the distribution of 
the magnitude of difference between present course and course 


from present position to destination, measured at specified 
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times during path generation. This measure is discussed in 
Chapter If and is illustrated in Figure ta. Recall now the 
procedure tor generating each leg of the tour. The direction 
of the randomizing velocity v) is selected by making a draw 
from a uniform (8,2%) distribution, and its magnitude is 
determined by multiplying the randomizing speed ‘determined 
from desired target speed through the water and drift, az in 
Figure 2) by the time Jength of the leq selected from = an 
exponential ©X. distribution. While the distribution of the 
direction of the randomizing welocity may be uniform, the 
distribution of the direction of the velocity through the 
water 1s not. The solid spokes of the left circle in Figure 
16a represent randomizing velocity and are equally spaced 
angularly. The dashed suvuokes of the right circle in Figure 
Da represent the water ve:ocity after the appropriate drift 
velocity is added, and are not equally spaced angularly. The 
angle co > that any water velocity vector Oo» forms with the 
X axis is a function of the angle (9,5? of the randomizing 
velocity, the drift speed, and the randomizing speed. In 


Figure Sb; 
= ve = 
+. <. a Sa 


where the randomizing velocity required to cause the desired 


mean water velocity, as derived in ¢?) and (18), 1s, 


4 





because the mean angle between rancomizinag velocity and 
course from present position to destination is %U degrees. 


The x and y components of the randomizing velocity are; 


Xx, = r cos9 vee sine 


G Gg a 


and the x and y components of the water velocity are; 
Xx, = x, + d ui > Yee 
Now, tang, = ¥ 7X 


tane, Oye / CX, + d) 


= CCC E/N) sine.> / CCM SD cose.) + di: 


Ss) 0 


= SING, ue [cose, + VV Od. 


And finally, 


= q 4 ‘i 
2 arctan tsine, 7 [coso, + VG Vi vn I Cle 


Notice that the redistribution of courses 2:65 not a function 
of lambda, but is strictly dependent on the ratio a 
and Knowl, o, as 4a 


by mathematical 


Even having the distribution af 85 


function of 98 obtaining the density of 98 


9’ 1 
analysis is very difficult, but can be circumvented somewhat 
satisfactorily by using Monte Carlo simulation. It is a 
simple matter to draw several hundred random numbers from a 


uniform (8,21 distribution and then tranform them with (127). 


This procedure was executed and the results appear in Table 1 
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Figure 5b: 


Representation of Course 
Angular Shift 








compared with the statistics obtained during path generation. 


Somé significant observations are: 


1. 


The mean magnitude of the difference between present 
course and course to destination obtained from path 
statistics, and the associated standard deviation, 
agree very closely with the statistics obtained for a 


by the simple simulation described above. This result 


Was anticipated because the expected position of al} 


experimental paths lies on the baseline even though 
very few actual positions fal] there. However, and 
very importantly, Figures 6a, 6, ¢ and 7a indicate 


that there is a strict maximum deviation angle for any 
te ratio. This is true only if the drift vector 
points directly from present position to destination 
all the time. I¢# position is always on the baseline 
(where the mean position for the entire process lies) 


then a maximum deviation will exist and will be arcsin 


ae >, #&a fact which Figure 7a makes convincing. 


d 
However, during travel the target deviates above and 
below the baseline, and though there 1s a maximum angle 
between the drift vector and water vector, to this 
angle must be added the depression or elevation angle 
between the horizontal and the destination. For 


example, if the target’s x-axis distance from 


destination is 1@ units, as is also the y—axis 


Sis) 
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distance, and the target is above the baseiine, then 


the angle between the horizontal line going through the 


target and the line going between tarset and 
destination is 45 degrees. Now, the next jeg out of 
that position might form the max imum angle 


Carcsin VOY pd between it and the dritt direction. 
Then, the pointing deviation which we are interested in 
is the sum of those two angles = 45 cegrees + 
arcsin (VM Nonetheless, the mean pointing 
deviation angles and their standard deviations are 
remarkably similar for the simulated paths and the 
simulation of (9); where one set of statistics i8 for 
the true distribution of courses and the other is for 
the conditional distribution given present position 
equal to mean target position. 

No significant difference is noted between the mean 
differences tor left unclipped, left clinped, and right 
paths. This 15S evidence that not oniy are left and 
right matching paths symmetrical in distribution, but 
more notably that clipping the left side by 
constraining it to match aright side path has no 
effect on the frequency with which the target points at 
its destination. 


Lambda, as expected, has no eftfect on course 


distribution; only the ratio V Gee does. 
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Figures é6a-c show the distributions of simulated couréze 
deviations for high, medium and low ratios. Notice that for 
a high ratio, greater than one in this context and signifying 
tight constraints, the distribution loads on the higher 
values, while for ratios less than one, signifying laose 
contraints, on the lower values. The break point occurs when 
the ratio is equal to one. While this loading, first on one 
end of the distribution range and then on the other, seems 
odd at first Blush, there is a reasonable explanation for it 
that Figure 7? helps to illustrate. For high ratios the 
tranformation of 5 to o, causes a high proportion of the arc 


length of the circle described by the maximum we to be 


subtended by the higher values of 9 Figure 7a illustrates 


1° 
this occurence and also provides a graphical representation 
of why maximum values of course deviation are low for high 
ratios. As the ratio gets larger, the drift vector in the 
positive x direction accounts for an increasing amount of the 
water velocity, causing the circle described by the maximum 
magnitude of the randomizing velocity to grow smaller. 
Notice also that when Wa = Ss (Figure 7b) the maximum course 
deviation is 9@ degrees and all feasible values are equally 
distributed because they subtend equal amounts of arc length. 
When the ratio is any greater than one, it is immediately 
feasible to have deviations as great as 186 degrees, the 


maximum possible. However, the higher deviation values 


suBbtend less arc length of the circle described by maximum 
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Figure 6a: 


Simulated Distribution of Course Deviations 
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Figare 6¢C: 


Simulated Distribution of Course Deviations 
V4/V,, = 0.26 
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than do the smaller values, as illustrated in Figure 7c. 
Recall shough that these distributions are not actual course 
deviatian distributions determined from generating actual 
paths, but are conditional course distributions given that 
present positon equals baseline position (which is expected 
position of target over all possible paths). However, and 
most aimpcrtantly, the mean and standard deviation of the 
actuai path course deviation distribution can be. determined 


using this simple inexpensive simulation once V, and ae are 


d 
specified. Thus, the target Knows how much it tends to point 
at its destination or away from its origin. 

As detaiied in Chapter MII, the third measure of 
eftectiveness 135 mean square radial distance between present 
pesition and laseline position. If the density of the course 
distribution after the application of drift C#(0 0] had been 
obtained, then ELR, “3 could be easily calculated for any time 
making the appropriate substitutions into ¢€6) and (5). 


Again, however, the derivation of the density of 9, 1S very 


1 
difficult and not necessary Because the process can be viewed 
as ao random tour carried out at the randomizing velocity, 
merely shifted right or left by drift velocity prorated for 
cumulative run time. At any time, the "origin" of this 
process, viewed as a driftless random tour with velocity Vas 
is the baseline position. This way, BERN 3 is calculated by 


substituting ek and t directly into (4). However, the 


resulting figures can only be compared validly to the 
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experimental statistics obtained for the unclipped left side; 
the requirement for left and right sides to match adds a 
further constraint tnhat is not captured by a process that is 
merely an unconstrained random tour appropriately shifted for 
agrift. Thus, the statistics for the matching paths must 
first be compared to each other in order to test for symmetry 
and then compared to the experimental statistics for the 
unclipped left side, which themselves have been compared to 
the figures obtained from ¢4). All these comparisons were 


made with the following significant results: 


1. Experimental statiztiz: +tor both the left and right 
matching paths support the a¢sertion that both sides 
are symmetrical. 

2. Experimental statistics for the unclipped left path 
mean square radia! distance agree very closely with the 
analytical fiqures obtained by substituting ae into 
€4). Experimental figures tend to be one to three 
percent higher, but any ditference that smal! 15 
acceptable as sampling error. Mean square radial 
distance grows with time, and for a given time ¢ 15 
less for ogreater lambda, aS expected. Table 2 
summarizes the statistics for the various time checks. 
‘Cleft’ denotes clipped left paths, ‘left’ denotes 
unclipped left paths, “actual” means obtained from path 


simulation, and “expect” means obtained from equation 
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3 rile. 


(Ht 


C4) by substitution of appropriate a The no} 
the table occur because, at lambda equai to one sor two, 
checks were not made for all times listed. 
Nonetheless, more than enough statistics were obtained 
to provide valid results. 

The ratios of experimental clipped/experimental 
unclipped mean square radial distance were calcul ated 
for the left paths and are summarized in Tatie 3. 
Linear least squares regression of the ratios aginst 
time indicates that mean square radial distance 
reduction caused by clipping (requiring a match with 
the right side path) is independent of V, and lambda 


d 


but dependent on time. The mean reduction is about twe 
per cent per time unit with a standard deviation of 
approximately 1.2 per cent over the range for whick the 
regression was done. On this range a linear f14 18s 
quite good, but notice that after fitty time units the 
linear reduction results in a mean square radial 
distance of zero, which 1s ridiculous because it is out 
of range for the regression. Cleariy, the reduction 
caused by requiring the right and left paths ta match 
must be calculated over the appropriate range. Though 
two per cent does not sound |)ike much of a loss, it 
becomes quite significant after awhile, as will be 


demonstrated when the random tour with drift compared 


with Brownian motion. Nonetheless, it is clear that 
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mean square radial distance for any path, given drift 
speed, randomizing speed, and lambda, 1s very 


predictable. 
B. BROWNIAN-DERIVED MOTION 


Again, the first measure of effectiveness addresses how 
representative paths produced by a method “look", and Figures 
Sa-g are representative of discrete Brownian-derived motion 
as executed by the second procedure delineated in Chapter IV. 
The paths in Figures 8a-d "look" acceptable; they show no 
pattern of regularity and seem to be executable, depending on 
the time and distance units selected. Figures Sa and 8b were 
gener ated from identical random numbers, a ratio, 
n® = 1.8, and maximum and minimum water speeds. The only 
difference between them is the physical parameter of which 15 
23.8 for Figure 8a and 18.8 for Figure Sb. The paths look 
Quite similar but do not look different enough for one to 


speculate about possible differences for the other’ two 


measures of effectiveness. It 18 worth noting that for 
a. = 5.86 there are more extentions and less clips than than 
for of = 180.9. This result 1s to be expected because a 


higher variance in the normal draws in €16) should produce 
successive positions that are farther apart than those 
produced with a jlesser variance. When mean square radial 
distance is discussed later, the amount of clipping and 


extending which takes place becomes significant, generally 
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causing tne ditfterence for that statistic for two different 
ot to be reduceci. 


Figures 8b and 8c were produced with similar inputs 
except fer lambcla, which is 1.8 and 3.8 respectively. The 
ditference between these paths is readily apparent, but again 
it 1s impossible from one example for each set of inputs to 
Graw any gereral conclusions about the effect of turning rate 
on other two measures of effectiveness. Figures S8d-g were 
generated trom identical inputs with the exception of Ye 
which is 9.12, @.23, 9.5, and 1.98 respectively. As the 
ratio becomes higher the paths straighten out as expected. 
The path in Figure 8q does not “lookK" acceptable, which is no 
surprise given ‘the strict contraints. Random tour paths with 
Similarly high ratios suffered the same straightening. 
Notice ats9, that as the ratio becomes high the process 15 
more constrained but less likely to obey the constraints. 
The gqgoai of the original controlled time of arrival problem 
was to assure that the target started from an origin and 
Visited the destination at the end of the stated time period. 
It i¢ clear in all four of Figures 8d-g that the constraints 
are violated, and furthermore, that as the Va p ratio 
increases the constraints are violated more. For any path 
that does not end at the destination there remains’ three 
options. First, it can be rejected outright. Secondly, it 
is possible that the requirement to reach the destination 


exactly is not as important as ending the path on time, in 
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which case, the path is acceptable if it is not too straight. 
Lastly, if the time constraint is not as important as 
ensuring that the target visits the exact destination, then 
more time can be given to the target with which Bie travel 
directiy to destination if the distance is not great. The 
problem of control will be the subject of further comment 
when mean square radial distance from baseline position 15 
examined. It suffices at this point to recognize that the 
problem exists and will probably affect the other two 
measures of effectiveness. 

The second measure of effectiveness quantifies how much a 
target points directly away from its origin and to its 
destination. For the random tour i: was found that the 
amount of pointing, as measured by the magnitude of the angle’ 
between the target’s present course and the course to the 
destination from present position. was a function only of the 
ratio nan and not a function of time or turning rate 
lambda. Those results provide a starting point for the 
investigation Of pointing for  Brownian-derived motion. 
Histograms of the deviation angle appear for high, medium, 
and low OY ratios for Brownian-derived motion in Figures 
Ya-c, and Table 4 lists the statistics for three o. values 
and $1x ratios. The table also lists values of pointing 
statistics for the random tour with drift obtained during 


path generation. Significant results are: 
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The amount of pointing for Brownian-derived motion 
agrees very closely with that for the random tour with 
Grift for low VV. ratios, as exhifited in Tabia 4. 
This result is expected because, regardless sf tne 
method chosen for generating paths, a drift velecity 
will bias course headings in the direction of the 
destination. Brownian-derived motion with the seme 
EON ratio as a random tour should exhibdit  simiiar 
Pointing behaviour. 

Table 4 also shows that as the drift/randomizing ratio 
increases the pointing statistics for the two difrerent 
processes begin to diverge. While the ransom tour 
statistics continue to agree ciosely with thase 
generated by the simulation of (€17)9, the Sraweriian- 
Beet et motion statistics decrease to a mean ot 
approximately 356 degrees and standard deviatics of 
about 26 degrees. Here, the effects of loss of control 
at high drift/randomizing ratios appear. The peinring 
statistics are more favorable, but at the cost of the 
target breaking constraints. Again, imagine the stick 
with a piece of string attached to each end as 
described in Chapter II. If one end of the string i158 
loosened and required only to be near the end of the 
stick, that is equivalent to adding more strinq and 
thereby easing the constraints. When the mean square 


radial distance for Brownian-derived motion is 
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examined, the lack of control at high drift/randomizing 
ratios will become evident. 

3. The pointing statistics show that for Brownian-derived 
motion lower values of o5 produce slightly more 
eointing. This phenomenon is most pronounced at high 
Va ratios, where the process is not under’ control 
anyway. However, it demonstrates that even at the same 
ratic, a smaller variance tor the position sampling 
causes the process to be more constrained. The slight 
Gitference in statistics would most likeiy be greater 
if low variance legs were not extended more and clipped 
1ess than high variance ones. 

4. The histocrams in Figures Ya-c show that the standard 
deviations of the distibutions tend to increase and 
Gecresse with the means. The low mean distributions 
are peakier than the high mean ones, indicating 
processes that are more constrained; the opposite is 


true of less constrained paths. 


The third and final measure of effectiveness for Brownian- 
derived motion is the mean square radial distance between present 
position anc baseline position, measured periodically throughout 
path qeneration so that a qraph or set of ordered pairs 
€t, ECR,.J> is produced for comparison with other path qeneration 


tz 


methods. Siqnificant results for this measure are: 


82 





ee 
BNP O Dat 
= 
— Lm 
ee 
= 
IL erat 


eos 
36 


Q 
us om AT 
Urdu 

meas 


<s 
Mg 
a 


L 


v 
NING O06 
ri 


Uy 

a 

LO. 

Wwe UM 


—_— 
foun 


Luss ly 
weCcrwot 
ot mt FP COM OY 
CO met 

usa Zo 

‘LE 

ota? 

Ju.u a 
er Be 


- ere ee + 


Wj teu aUtuas 
ceeds a a 


Ce en ee 
ee @#e¢ee 
UN ALS mt N 
~ LI O 

eo ee @ 6 
US AU UY 5 arin 


a aca 


ONG mt at tm et CY 
QONoOCcIdSGdy 
Oeet Vt OD 
TarOoogny 
POLIO wiser 
mets OT 
CICONLAE NI 
=e PiNG Tas wt? 
er tire + 
LL tt uy 
PP IODyD 
S&XLursctaM 
INCRE ONT 
NN me mm OS mt 
AGcDeaos 
Se ie 
AMNATEAH MY 
NMmOC MU 
Warctigmme oO 
Tt (CMthate 
— 

Men tet 
MOlLCE AD 
-~aeHeoeerhe eo 
AL SLL atts 
ON Ouratwo 
ACNDINGD 
™ OO mt QP AIO 

- Fampv~yv Ow 
TUVATOrWS 
Veer pwo 
Oe OW Cw 
Loko rm OL et 
Ey Sr 


wiki 


WOR INN TW 
NOILMAIMYISTA 


—e 
< 


Mw ToONMOor vo 
re OACOCO 


3 
j 
4, 
» 
vas 
€ 
I 


cm ao 
( 
Ww =>wM 


2» 


3) Sao 


ial 
uu 
= 
GC 
oe 
3c 
oO 
tJ 
rv 
qa 
= 
us 
|— 
ab 
<t 
o 

b~ 
ih 
% 


"97 "89 "e "19 "94 49 “¥4 a ay ft "27 an *y 


CHER EEF ETE TE Re AER GEER E HER ERE CHE ERE FE REET HR RHE ENE NE KEH Hee SHOE H He Ke meee CHS ee eet ER RE tee 
SIAGESE ALR HER Bee KES MEE Ee HK TH Hot CHK HRE ENE HER Kak OER BER te FTE RvR FRE See fe 
SHLAA Rm Bee HEH HAE HEE CEH Aad HOR EKG SHE ONE Ree Oe oR AEE AHH ThE Eee LEE EHH HHH 
Sim wet RHE HAR SAK HEF KR Hee HEH 

SAY G AXE HER HEY HEM HEE HEL HE 

SASed4 wm4d KOR HEH Hak HVE 

3 24-2 stk Fe 

a eek FON 
FVE 49% 

4 mee 06 

eee Nk 

Ate ane 

for ee 

eter tjgee 

3 ate 

wae a Se 


me me al vale . 


ite “333 ee 


e 


* 
* 


ted 


take Ue Mev JG 
RAK HHE FES 
Swe 2% ARG 
eat @€%5 


Gee Na eK Tee Pete ahr aes ; 
Y 
4% 
é € 
* *4443¢745454 
x 
4 


Ps 
OR Wel eee ek O Roe Re RR F 
XE CNR He AEE TC HRaR C1 t FEe 
Sen Gale re) epee Ree 2&* Isom Y 
1490 bay wth eer eer *o¥ toe + 

a 
¥ 
€ 
4 


ob at 


Ree BEL THIS Hee why 4 3 
Sue Fee aki eee tke net 


eee KER RHR KEE HORE KET 


emt F He 


able kth €h eo vee wEi ste 


eed 
MES 
e?A 
rae 


to 
prim? hum mye’ vr & deg 
ES SOe tT Ee ek 
oe ere? eke e bf 
: Sky BE KN AH Oh 
far ath Boe OR KER RED HUE 
, 
¢ 
t 
& 


wg we we ee we & 


Sete GS Tae 


[O24 ght?" €t€ +Ot 


| 
raves 
ane 4 
ea 4 4 j 
eile? Oe wee? €4e @ 
CARTE ORE DADE *& 4 
thia See #€Oo° et wad 

SOF HL TSOP REE 

ek ny e° tan €£E* 

+ 


*i. 
eevee 

ett Aw tn eet eH 
Pw +e 


wt te wee Rewer PB ce ce 
ww & «rp 


**¢ 


i] 8 ig 


Lee Lte eee eee eee de et +H 
wu ef Heo x oF f- 4 
ie + 

Le @&€ tH PH See mae eE HF HOS e Dy H 
eenxenwxta eter Pee en wet Ae pe ee t 


lee 
te 


& Putte anes 


ag ° 


we Set ot + 
Ket ay £e#ewea & 
od 
= 
So 
° 
<- 
> 
> 
bd 
te 


0°OL =,0 
. 3 18 
3) ae 


UOT}JOW PeATAaep-ueTUMOUdgG UOJ 
SUOT}eETAdG asdnog jo uotynqtd sta 


Ol 


>PG suNsTy Z1.* 


eeeeerte#eeeteeeeetkte#ee?#e?ehe @.% © @ @ @ @ 6 


\! ra T 0 e+e Z 9 6 Gt 61° 22 12 6 62 19. GH °SCTIOTIT 1S1 O88 ole [02 27EI FIT 9£ 95 Sete Gees 
OSAT = 371S JIWdWVS SITINSNO S34 


8-3 





Commotion — | iiUuphe Vy TQe*yguvHogbai ed NV abe bd be 
i RA te  S 11 LNG 2” 3048 F195; DUE c€ ~Yia TOV rh I4G PO | — IVA NWS OL Vee, MY Gn kilav 
Oey betecr il CHOU JPULINVAO see C9 VS AV *f (Vidgu 7 (yeas | JORN TO + _VOCLOE Od av eVUUE 
We gsoneeS°S CNVIGSW) JVILNTAS Q9° Oe gh eulZu"y S1soLren» 10+ 376-21." c AIG WV aw 1044526195 °S Hv ot Glow 
TO¢+¢Agec 153°C (JIONITH) JVLINVIO Se Q- ss Fh ok EERE (N-—FIGOG OLE NS HVA 43uUz TU45/79C594°S NNVaeloa 
1O*3810225°1 JFVILNTAYD ore ON SOPPGe yore Lal | 104#4GUGZO°E AIG als TO4 5eGGLtS°5 NI1dsan 
SOKAVSPVE ITE WHOWINT A O44 99076 F6° rae} 2O*sei sg? °o JINVINVA TC4+374%6C9L °S NYiW 

NOLTLAGIMLSIO SINSkRIa IWYINSD AWIHNOTH Q0sudd5 ADNIINSL Wuil'tasd 


“JLT ao ee Br i i ee oe “83 oa if eal Be a a, we aT are 
AERO COC CAST + CHEERS GOSLSTEEET EMANEG CIO 4EtOA Lt REEUE K} Fame EOE yAl cEhea OKO” SPAT EVEL SPN EOE SE eee tae FeO Sth EO Ee EN eet ee rey 
soar Sp4d4eme Jot dvee C4 SHE HEH KE HK Kk COR HEE ROK KEE Hk HI HOR KEE HOH til® Oe KER Ho CHE CHO HE £O4 CEE TER 4, 

eh, S44 S3HG 1 AAGZAAGEE BER HEE HHH KHER KEE LEE EEE HEE KEK TR HER HEE OHO BEE He OTS RMD HHO ED KEK OWE HEE Bde 
**3 4 S44 BER HEL VER HHH EEE REE Kee RE EH OR OSE KEK EW EE ERE Heo Vee Ree eH RHE HE Rew KAYE 
Sd KOR EER KER HEE ORE FEE NEE KEE OTH HOE POE OH” Bee ed OT HE Kae HH EK TER CFE 42 
44 weIS4dS HOw BH KOR Ee KER KE HOH KEE EHH EW? TEE BEE HET ERE HER EER Bek PLE wee SEL TO® 
Jide A eeK EEK EK REE HT FE KOR REE Hem Wo Oe REE HHS RK bee AO FOr eee REY aes 
S¥SACEK ERE CEO BOK Hk 7H EOE CHK EN FER Ree EHP KET AEE EO FER EK HHH HE 
S34 kh HOE EER Ree CE OH THe BH CAS BOW BKK HET OHb LER EAE FEF Bee XE HHS 
KAwv whe HEE AER whe woe RHE EOE HN ome EEE HET CHE RH HER HHH HHH SEE FO 
4d Ee Bek HTK HOO BH HRT ENT HER ERR HHT REE FEE EV Fae HEE BE FTE 
dd ate bik €E PR FAO EHH ENT FO FEE KE” COE NER Ome Hee KEE BEE EEsICV” 
KE O68 9G £°R OOH HE CH” CHE BEE Ro HEY 444 CHO ELE CER CHE HO? 
4O7k Apt 4 wld EAR ShE EW? BEE HEE Fh HEF TEKH HEX HER Het EA 5K WHA 
ke BER OE Om HE ERE Oo Dae wee HES Hee KOE ES FOR AEE DAGdte 
SIS CHR HE POH EME ERO EW? HOr Hie vem Hee BEE HHH et HUE 
WED AEH EEE HTH CHE HER EN HED Few KT Hee ANH Bee akT FEED eo" 
SEGRE COR ECR CHE EE EWS He FRR Ro HEY HEE Hee Hee 4A 
SSE KER OE CHE THRE CHS GOH HE O° EK DEE EEE BTR FG 
Sem FUG OR CEH FOE Chl? VE HEE De? WHE Ke whe BOSS 
eG EEE 7H SEE BEE OAT HER HEE SH” Hee Bee 4iwm FIE 
STE RE ETL ORD KEK Ol” FPR wkE OF” ERO KeHEGeEG J 
eR mee TOE EN VET EN? Pee KEE EHS FEO VHT 3 70° 
3S HE wR wok wee HH Ree ae Em ER OF 
weer °R HVE SEH Ne wee EER Ho KOK 
Sux w°R oe FEE EN? ECHR CKE EHS HEH J 
44 3d OPR HOE HEH BH Fey COR Oe bees 
eos Per tre EW Fee eee 2° CoM GO ° 
4 ow wey Hee eel? OE Oe EF BAe 
€°R tor Few KWo wee HEF te? 3 
e°R FOR BH ENS CRE Oe x? 
Ae°R Oe CHE ENS eK Here 2e® FG 
m°@ MER COR FT Ke HEE GO? Iu ° 
3° €¢a COE OH” Vee ee 7° 
c Cds tee We we aye oh tae? ZF 
O-O eae) 4 Jae HEH EH KKH EK Gm? 
a od i J 4k EN” EK OE £¥® 
° A4eee On cer woe tod 
I Pp : dee oW° r41e €ee #5S LA ae 
O°T - A/ A Re CW HED SKK KD. 
oe y vk OV Ow RE . 
: ww° @er ee Jd ° 
e e & e 
UOTLOP! PeATdep-ueTUMOUg TOJ + aera > po" 
7 eu? Sik tee Z 
SNORE TACT 2SihOp 97O UCTING rays tq . 3 wed wee 8 
‘ Ad Jueves 6° 
‘ 4jd° wee ° 
-qG. ound t J : Hr dee bo * 
e e ° € e 


t iE S 8 E Tr SG OT Tr et 22 T2 O€ 6E 22 99 46H EOT BET 957 69T GPT ZaT SOZ 5OT BIT 16 28 48 89 9G 65 
0622 = 3718 sI3WVS SIL InINGS HS 


84 





209975955 ° ft WOWIX Va OU4A3OeU6T°S NV sn navn 
CO+FIF29719° T JTELNVAD 05° ros eda a en ei eat OVANGSOLS. ye eee NVoh bY 
€04326959T° I (Z9ONIH) 4TILNYNO S2° 904#551G50C6°E TVids 70+3L119652 °T 2UNT 4 #3dC€POS°R ION Teja 
TO#9E5Z909°2 ENVIG3SW) JIELNYONO OS? TO-29uUSt 80 °B- SISOLYINY 1043602506 °c AQU NVan TO *21c656°2 Toh be 
TO43D9Z0ES°E (JONTH) STILNVNO SzZ* TO-sc66USL*E SSSHMSNS TO-3SCEOCT°9 VA 43d) TC +#a1SP 552° WV ah bao 
TO0*365E5e27°T JTIANTAD OT? 2O4sTcotzo0°t UW TO+t3BTLZ20L°4 AO GELS 1C+3C6<¢9CO0°L NhylUvdn 
cO-35725TL9°? WOWINIW 6O+3Tsl2G0’°e | cO+37SSTIZ°2 JDL YUVA TO4529Z2TL9°%R NYAr 
NOTLASTSLSI90 SINSWOW FVHLNID J3ZHOTH Ovauds AINIJINAL WolNse 
"O21 at a i: tf Mey ace a a Oy es “a 7 ay ey ‘Si a oe 
MESH HEREEEEFSERFES EF EEF HEHE REESE E REE HERE ERD RESO ERIE MEER EER REDE REN OES HER ERE ERS EERE SET EK He EEE RTH RE SESE SAREE CES TIORS 
HEE REE CHE HEH HEH EE FOK TER ERE COR HER EKE THRE BSH BER HEH BIE EEE Wee OOM HER HR UKE TOE HET HH Hed She WES CHE HOH Fr 
EOS EE HEE PEK BEE KEE CHK KEE KEE HAR Hee KER CHE Eee EE ERE Ht HEH WEE wos het FeR Et EEE Ho FOR OSE KAR WK HEH HER FTE 
SER CER KEE COR EE HER EER CFE HED FOE HEHE OH CRE HE OEY FE PEK VER HED KOR HEH CHR HET HEE O° COE 868 BHE HEE ChE RE tOH 
SES HEE KEK KHE CHE KER HEE HHA BOHE HHH BEE CER SOM Hee FEE HEE BSH HE Nee Boe WKH FE wee EF UOT ee ETE EH HAH ae BKK 1 
eee THE HHO ERE SEE SHE HEE EKG FER FEE HEE KOR THE EHH KO AE Bek CHE Wee K°R CHE CER KER ERE KET FEO VER HER HIKE HEE Hek UL 
EHS KER EEE EEE HHH SHH HEE CHR ERE HEE KEK SHH CRE HEH THE SHE AEE EO WEE °C HEE KE HEE BEE HH” CET VER EHH AK KeR KER HL) 
CHE KER HER CHE HEH Che KE HRHe KET HHH SHE HHH CRE HEE SEE HEE FO HEE Wee FOR HOH HOR HK BER HES HEE oem Cae ERY KEH KTF Ses 
HE FER FER EEE KEK CHR ERK CEE KEK HEE BEEK ERE ORK Et Hes Khe KR OTE HEe 8° ENE CHE OOH THE HES CHE CHE HOt SEH OF KER ZTE 
tee KEK HER HER OEE RK KER HEE HEE EEO 4k FEE ORE SEK HAG HE REO HEE HOH FOR Bek COR HE KEE KES wee HER Hee FET HOR ERE Lok 
KR EH OH SHE CHE BEE CHE HK ECE CRE FEE Cee CRE SPH EWR SES BIR BALK Het OR HEE ERE HER €OT HHo GRE RHE Hae ete pew awh see} IS ® 
CHO OOF COS CHE HOS OK FER KEE ERS HER DHE BOK ORE HHH CFE Fee WE HhE Wek 8°R ChE Oke HEE CEE HET HED HEH HET REE KHE Bee HOR 
Feet EC OK See CHE Sak FEHR FEE KEK COO BEL Bae CHR EEK KET Ao WHEE COR Wee ¥°R Oe OHH Hak ORR Fo EE DEE a2 e PEA HEE bet Wee 
EEE HEE EHIME E HOO REE COR HAS KER Sek OEE Hee COR HA HEE HED BHE HHA WHR OW HEE GET HER REO FES EKA ERR HK ENE CER FES wre 
2G ChE HAE HSH FEE It HER CHE REE HE KEP BEX CHE HED KH UKE EER HEE Wee 8° GOR LEK BHE Sew b> BRE SHH EMH FIR THE HRHE ae 
Oe 4H OH Oe HEH OHS EE CHE HEE OEE ERE HER CHE HEE bet he She HEE Wek £98 EER EEK OME KK ETO bee AER FEE HEE HE tHe veh 
4 EG OH) AE CH XO ERK CHR EET HAH CHK OEE COR HER Hee Bee EEE Ree Wek oe HEE HEE EHF FES HHS KEE TEE HER HTK bE HER BEF 
d4% *e43544 SRHASEHISEK KEE HYG KEE HEM HEE KEK OHH CF4 ERE EHE CER CHK NE O° R Ae CEH CFF BEE HHT HEX he wee FOF BH wy > Ot 
ethijdee wee 4 Skee OR EUR BERK HOR tk SMP Fad KE FET RTE Hee Wek 8°R Beh He SEE KK be Vee Nye KE+ HEL HeR Fa FUK 
ae 4d jad Se EE EH CHE HEE EE “HR BHF SHE EE E48 £4 Wer HoH BE HHH THEE HOH EHR” Het HEH SLE KEE CHE Beh tet 
SWERSSSAA AR eSIARE SKE Kew CHR HOR HE FED EE KEK HOEK FoR NH GER KEK KEE KET eke FET EEK HTH HET HEV Bred 
J Wee SHE REE ESE HEE Hee CR Ohm HEF Kee FER elk WHE O°* Hee KEW HE HE Reo be Lek wee yh ee EH BEPC ® 
SAS OH HE EES HE SHHI° SS CEH HTS HEE TE ASR Wet HOR BEC HOR Bete KER EHO Ket CRE TARE week EE BEE Es 
See OH SHE ER Hw HER Het FRE €8% Bek Wet FTE BEE KER HER BEE Te” Oe Che Ee HHA ERE TEE FO 
4 J ER RHE ORR HER HkHE HEF REO KER Was °R TEE KER FOR She BH” SEH ERE HHH HEE BER BEE CAE 
S3¢jseed + week KOK Cee FH EKE HEE OR HEW EEN FER BRE ROT REX wee ve Eat Lee new 44? 
3 SOJm OER THE HEF EEP HEE HHH KEK Woe ao KRER HOKE HHK HER HET BKR HK TE Ee Het BEE OE? 
J He "HH HES SOR Cam OK RTE WH °F CHE HEH HD Hak OT HEHE FEHK FEE HHH HEH BRE Fld 
SHO ASK JER ERE HER HE BHR Wee FOO ERE HER GO Hehe em CEE EK HR ThE THe Har Hike 
S$ 7 Se SIG RHE FHE KOO OEE Wee or KOO Hee Ok BEE HET Cea ERR HSK KEM HEE EO O44 
‘ Pees eee IK HER ETE HAE He KO HER COM HEED CHK FET Fee BEF AF JAGHE EHH HOD a | 
a. Sd Vee He RHE WEE HOH EER BRE REA AEG CBo Ld CEE 4H HIE Hoe EHEAAde REE? 
°J4 Rae deee Se HOE WER BOK HE HER FER CHE HEo KEE FOE SAIAGHHE FHL BEX 
: 44 Jed bed WEE EON HOH HER BHEx hee FES KEE tt ad ev dik te 
0 "OTL = = eesdese sek Wee 8° O KEM Hee Che HEE CHF VEE the eee BOR 
7 D . eS Htedede Nae EPR Eee OAH HER Oh OH” BESH HES Bed he 
x4 a %4I% $24 Vite OS Ree Bek Gee ESE ERT TTR EEE Rte ¢4e 
ee etedietde’s bee #69 BED BEE HET CN O4H ee 
; A Pp ‘ 3 ee Wie 3° s wth bee SEH HOt FE bev SHH 
SC Ose heen : Wed &°S FER Ee OES Ore eho OTH Rae J 
i; W 6° ¢ Ie HEE HH HED BHT KEE BH 3 
i Ww woe ED He EEK SERASE® 2d tee oo.” 
UOTTORP PCA Etop =e TUMONG., AO W Settee eee CO EEE HET SEF Hed 
° ° ° W + is 4+F¢ a ete tes #0” sHe 
W ei Ae BEM «HE RE KEES 
uoTIeTASG SsaNnoD jo suoT INndtTaisTd yt) bah tees Fem aceibey tes 
. y 44 y- rhe ‘ 
: ee 24@ vee i 
[0G edhsTt yj W ° deasdize ade 
g Ww ee 4a% J id 
Re mm trem Fem Fae ee teen be me tren Fem eho me hmm ambien oper eterna teen ta wn town tend en Bn tn ta a mn Ba a tae nm tan dd we a ben net mr eter fpeentqenrnt--=+ 


Rel@ocdt OF GUT Veletct ESI <SIl clLI-6St 26 102 Ele Ole Les 192 Osc B9C 9lic FOU GIE THe Sve Ger WIC Tet fhe Ocd Ste eae C 9c eee 
OSZL = 3218S JIWWYS SSIIN2NO5a4 


85 





paw 
a 


Figures 14a-c illustrate mean square radial distance 
for four Ned ee ratios and XX = 1.03; each figure has 
eurves for three different values of the physical 


parameter 6 Immediately 1t 1s clear that the mean 


2° 
square radial distance decreases as the 
drift/rancomizing ratio increases. This result follows 
directiy from the ogreater constraint placed on the 
target at higher ratios, and is similar to the results 
for the random tour with drift. 

The problem of control for Brownian-derived motion is 
apparent in these figures. Even in Figures 18a and 
1@b, at the lowest two drift/randomizing ratios, the 
mean square radial distance does not go down to zero, 
indicating that on the average the target does not 
visit the destination. Notice that as the ratio becomes 
progressively higher in Figures 16c and 18d that the 
mean square radial distance tends to become even 
greater at the end of the time period. Clearly, as the 
croblem becomes more constrained, control is lost. 
Anather characteristic of this path generating method 
with small lambda 1s the hook in the curve at the end 
of the time period, indicating an undesireable increase 
in mean square radial distance during the final time 
unit of travel. This increase is the result of the 
truncation of the final leg when time runs out, though 


the random numbers for the final leq were generated as 
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if no truncation takes place. The simple soluticn 153 
for the target to travel directiy cowar the 
destination on the last leq in order to get as ciose as 
possible. A target would probably do that anyway since 
1t does not really need to be told what to do on the 
final leq. The algorithm for path generation was 
purposely not modified in order to demonstrate this 
peculiar phenomenon. However, it is importan: to 
realize that a true “bridge*®, Brownian or otherwise, 
does not exist as long as the two endpoints are not 
connected by the target path. 

At lower drift/randomizing ratios, the mean square 
radial distance is less for smaller Fos as Figure iGa 
clearly illustrates. However, as the ratio inereszes, 
the mean square radial distance 135 approximately the 
same for all three values of o5. The curves in Figures 
16éb-d also show that for different o., values the 
maximum mean square radial distance occurs at dit¢erent 
times for a given drift/randomizing ratio. The mean 
square radial distances for different on became closer 
in value because of the clipping and entending process 
that occurs to Keep the target within speed limits. As 
the problem becomes more constrained at higher drift/ 
randomizing ratios, the target tends to bump into the 
speed bounds more. Hence, the process is clipped and 


extended more and the natural differences among the o. 


Pi 





values dominate less than the speed limits, which aii 
paths must obey equally. Also, as the mean square 
radial distances for various siqma square all get 
closer to the baseline for whatever reason, they are 
bounded on one side by the baseline (the minimum value) 
and get sandwiched together, as Figure 19d clearly 
demonstrates. 

Figures 18a-d show that for a given Baa there 1S a 
value of sigma square which produces a mean square 
radial distance curve that is symmetrical about the 
line x = T/2. For Figures {@a-c the values of sigma 
square which come closest to symmetry are 16, 5, and 1 
respectively. Notice also, that in Figures 14a and 
1@6b, in which the drift/randomizing ratios are low and 
s$iQoma square equa! to ten and five are the values 
closest to producing symmetry, the curve for sigma 
square equal to one 1s very asymmetrical and tlattens 
out for big t. This occurs because for these low 
drift/randomizing ratios, o. = 1 over-constrains the 
process; it forces the target back to the baseline too 
soon. 

Figures tia-d are different from Fiqures !@a-d only in 
their lambda value which 1s 3.6 instead of 1.4. Tt is 
immediately clear that an increase in lambda causes a 
marked decrease in mean square radial distance. A 


Similar result was found for the random tour, for which 


ag 





an increase in lambda caused the distribution of 
radial distance to start piling up around the expected 
Positian o- the process. There are some other subtle 
differences also. At low ie ratios, all curves for 
the three sigma square values are skewed slightly less 
left than for higher lambda, and for higher 
Grift/randomizing ratios they are skewed slightly more 
right. Notice also that the curves for the lowest 
drift/randomizing ratio in Figure titla are closer 
together than they were for the same ratio at 
mw = 1.4% tn Figure 1{@a. This, again, 1s the result of 
aljl curves being bound on one side by the baseline, and 
therefore forced together more at the lower mean square 
radiai distances+e which higher lambda produces. 
Fimaliv, the "hook" on the end of the curves is gone, 
most likely because the mean number of turns during the 
final time unit is three instead of one. Consequently, 
the target turns a few more times during the last time 
unit ot travel toward the destination and thus does not 
rely only on one leg to hit or miss. Again, however, 
it is likely that a target would travel directly to 
destination at this point anyway, unless it desired to 


adhere to its course change policy strictly. 





SWIL 
o°s9 O° oss O0°0S O'Sh O'Ob O'SE ODE O's2 





o'oog =sotost §=sooot.sS sisi 
SONULSIO 


0°0S¢e 


Socee (HH 
~---0'S 
aaa Olas 

Gee 


0° 00% 


O° OSE 


plo 0 


0°O0% 


NOITLISOd SNIVSSbE WOYS 
JONUISTO WWIGYY SYYNOS NYSW SPT eansty 


94 





O° OL =e ees eee: 


aC 


0 
0° 
0 


Wid 
Ose Sez WetO2 S"2t ost s zt aor 9°z 
3 ES) eae) oe RE ee Is Eee ees Oe 





° 
Neveapnils o_o. nciWelets 
ee 
@egeet? 


NOTLISOd SNIVSSYS WONS 
JONBISTO WIGKY JYBNOS NUSW 





0°0S2 0°002 o°0st 0°00! 0°0S Oh 
Biiethsy (0! 


0°00% 


‘(TT ean3sty 


ws 





o°s! O°sT 


WIL 
o°21 O'O! 0°8 0°9 O°> 


85 °0 


NOTLISOd ANIVWSYE WOYS 
SJUNUISTO WIGUY SYYNOS NUSW 


Q°o 


\ 
o°szt. a°ost: «=arsz@ts OGG sssaseC OS Ssé'P 
SJINULISTO 


0°00¢ 


PO ocd & 7 


es 





0°6 


o°s 


ed 
ao®* 


SWIL 


aeeecne Vr enmremrenwaan eo 


| 
t 
i, 
t e 
i 
2) 
=) ° 
>> 
“i et 
| | 


SS) © 
om FS 
it Ut 


-~ 
7] 
{Cc 

e 
oN 
{! 


NOTLTSOd SNITSS8E WON 
SINUISTO WIGkY BSYYNOS NUSW 


0°O 


v°O 


o’sét 6 gras! «=oszt.—s cok Sask SS Ss 
SONYLSIG 


0° 002 


we LL odie TF 


Sa) 








VI. CONCLUSIONS 


Chapter I! defined a controlied time of arrival scenario 
in which a time constraint is placed on a target that i5 
required to travel between two endpoints. It 1s desireable 
for the target to "randomize® its motion during the transit 
in order to provide the enemy with as little intormation as 
possible about target origin and destination and to make 
target detection and redetection difficult. Chapter I] 
discussed desireable qualities for such target travel = and 
delineated three measures of effectiveness. against which to 
measure any procedure for oproducir.¢q¢ terget paths, and 
Chapter III examined previous investiaqatiors into the problem. 
which have provided direction for the two approaches adopted 
in this thesis. Chapter IV deecribed the random tour with 
Grift and discrete Brownian-dervied motion in detail, and 
while the results for each metned were presented in 
Chapter V, the performance of a path producing procedure 
against the measures of effectiveness 15 not important in 
itself; but rather, the measures provide away to compare two 
or more procedures, one of which can be judged best for a 
Particular situation. 

The amount a target points away from orign or to 
destination is very nearly the same for both the random tour 


with drift and discrete Srownian-derived motion. Tables 1 
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and 4, which list the mean values and standare deviations of 
the course deviation distributions, show neariy identical 
values for both procedures at 1ow be ratios. Differences 
Begin to show only as the ratio gets bigger. But thouch the 
statistics become more favorable for Brownian-derived motion, 
it 1s probably because that process is not under comosiete 
control of the constraints when the drift/randomizing ratio 
is high. If 1t were completely controlled or if the ransdem: 
tour were also allowed to violate the final visitatiaon 
constraint, then the statistics would probably be very near?y 
identical over the whole range of Y ratios. 

The failure of Brownian-derived motion to contra! targqer 
travel by quaranteeing visitation on time at the destination 
‘15 a potential weakness of the process; it literally does nat 
do exactly what 1t 1s supposed to do. However, the weariness 
is not important just as long as it 153 not vital tor the 
target to get exactly to destination exactly on time. OS a hs 
is vital, then clearly the choice is in favor of the random 
tour, and none of the measures of effectiveness are retievant, 
except for the requirement that the path “loak good" and be 
executable as described in Chapter Il. 

If the target does not need to have a perfectly 
controlled time of arrival, and can either fall! short of its 
destination at the appointed time or else take extra time to 
get there, then the mean square radial distance becomes’ the 


important measure of effectiveness; pointing to origin and 
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Gestination is similar and representative paths of both 
procedures “look good" and are executable. Figures i2a-d 
lifjiustrate the mean square radial distance between present 
and baseline positions for two sigma square values of 
Brownian-derived motion and for the random tour with drift 
with and without the requirement for the left and right paths 
to match. While the mean square radial distances for random 
tous paths which do not match have no practical siqnificance 
Since such paths meet none of the constraints, they do show 
how much radial distance is lost by requiring that the left 
and right paths meet. Recall that the loss is about two per 
cent per time unit, and while that sounds low, it is DSecomes 
very significant as time progresses. The fiqure of two 
percent was obtained, -as previously described, by Jinear 
rearession and is valid only over the range of time for which 
the regression was done. This limitation is illustrated in 
Figures il2a and izb. Notice that the curve for matching 
random tour paths begins to saq in the middle, when actually 
the mean square radial distance should be increasing, however 
sliq@htly. The sagging is a direct result of the linear 
regression operating at the edge of its valid ranqe. The 
curve should be rather flatly rounded in the middle instead 
of saqging. However, the curves presented are accurate over 
their range, with the exception of the slight saqgqing as 
described. In Figures 1i2c and i2d, where the time period is 


shorter, the curves do not sag at all. Most importantly, all 
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four figures show how much mean square radial distance 15 
lost by requiring the paths to match. Under the conditions 
prevailing in Figure i2a, the unconstrained random tour with 
drift exhibits a greater mean square radial distance’ than 
Brownian-derived moticn does tor either sigma square value. 
However, the mean square radial distances for the matching 
random tour paths are much jess than those for Brownian- 
dervived motion in ail figures. Notice also that while = an 
increase in the ae ratio or lamoda each cause a reduction 
in mean square radial distance far both procedures, that 
increasing lambda affects the random tour process much more 
adversely than Brownian-derived motion. This is a point in 
favor of the Brownian-derived srocess. 

One might arque that 1% is natural for SBrownian-derived 
motion to exhibit higher mean square radial distances, if 
only because the process is net under control. This arqument 
is not compelling and scne need only to look at the figures to 
see the great disparity. It 15 not 11iKely that the failure 
of the mean square radial distance for Brownian motion to go 
all the way to zero at the destination is the reason that it 
is more than double that of the random tour in the midrange 
for three of the four cases. However, one might §arque 
further that the constraint which requires the left and right 
random tour paths to meet up causes such severe deqradation, 
and that Brownian motion might suffer similarly if it could 


be made to meet the constraints exactly. Nonetheless, the 
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paths perform differently as they are. and clearly the mean 
square radial distances are much more favorable for Grownian- 
derived motion if the strict breaking of the final visitation 
constraint can be tolerated, and remedied. 

Recall that the first measure of effectiveness which is 
applied to a path producing procedure checks representative 
paths to see whether they “look good’® and are executatle by a 
candidate target. Failure against this measure autcmatically 
disqualifies a path from further consideratian. While both 
the random tour with drift and Brownian-derived motion passed 
this vital first test, there is a subtlety in the way hoth 
procedures are executed that makes Srownian-derivec motion 
more desireable. In order to obtain matching left ard right 
paths, the random tour must be executed in 1%5 entirety 
before the target begins a journey. The tarqet must ‘chen 
follow the instructions closely to make alli the correct 
courses and turns. While this procedure is possibie to carry 
out, it is quite exacting. On the other hvnand, Brownian 
derived motion can be executed one leg at 2a time because 
after each leq the controlled time of arrival pronliem 15 
reframed as a totally new one using present position as the 
new origin. Hence, perfect navigation is not as critical as 
it 1s for the random tour. In a sense, starting the process 
over after each leg always gives the target another chance, 
Just s0 tong as a flagrant violation does not occur which 


causes the target to be faced with an impossible transit at 
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the beginning of some intermediate leg. This feature of 
Brownian-derived motion should make it more saleable to 
target captains. 

Both methods for producing paths should be available for 
a target to choose. Clearly, if the final visitation 
constraint absolutely must be met on time, then the random 
tour with drift method provides the only guarantee. But, if 
the target needs only to ensure that it arrives in the 
immediate vicinity of the destination or may arrive at the 
actual destination slightly early or late, then discrete 
Brownian-derived motion performs more favorably against the 
selected measures of effectiveness and is easier for a tarqe?* 
to execute. 

The two very different methods presented here for 
generating paths represent only two among many, and the 
variations on these two procedures alones are infinite. For 
instance, the random tour with drift could be modified so 
that 1t was executed from one end only, instead of from Both 
ends as done in this thesis. After each leg of travel the 
problem could be reframed as a totally new one, in much the 
same manner that Brownian-derived motion was restarted after 
each leg here. It 1s also quite possible, and desireable, to 
devise a way to force Brownian-derived motion to arrive at 
the Gestination exactly on time, satisfying al] the 


constraints strictiy. Thus, the two methods presented here 
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ar2 not only quite workable, but also provide suggestions 


for turther investigation. 


168 





1@. 


LIST OF REFERENCES 


Washburn, Alan, "Probability Density of a Moving 
Particle,” Operations Reasearch, v. 17, no. 5, pp. 86i- 
871, September-October i969. 


Daniel H. Wagner, Associates, Interim Memorandum to 
Applied Physics Laboratory/Johns Hopkins University, The 


Ornstein-Uhlenbeck Displacement Process as a Model of 
Target Motion, by Barry BelKin, Feb. 1, 1978. 


Daniel H. Wagner, Associates, Interim Memorandum to 
Applied Physics Laboratory/Johns Hopkins University, the 


Random Tour Target Mode! with an Arbitrary Course Change 
Distribution, by Barry BelKin, July 27, 197%. 


EPL Analysis, Random Waik Statistics, by Edward P. 
Loane, Nov. 4, 1977. 


Daniel H. Wagner, Associates, Interim Memorandum to the 
Applied Physics Laboratory/Johns Hopkins University, RMS 
Velocity Behaviour of a Constrained [QU Process, by 
Barry BelkKin, May 5, 1981. 


Daniel H. Wagner, Associates, Interim Memorandum to 
Applied Physics Laboratory/’Johnis Hopkins University, The 


Random four Target Model with an Arbitrary Course Change 
Distribution, by Barry Belin, p. 11, July 27, 1979. 


Daniel H. Wagner, Associates, Interim Memorandum to the 
Applied Physics Laboratery/Johns Hopkins University, RMS 


Velocity Behaviour of a Constrained [0U Process, by 
Barry BelKin, p. 8, May 5, 1981. 


Ibid, pp. 9-13. 


Daniel H. Wagner, Associates, Interim Memorandum to the 
Applied Physics Laboratory/Johns Hopkins University, A 


Result Concerning Datum Aqing, by Barry Belkin, Dec. 22, 
97a. 


Daniel H. Wagner, Associates, Interim Memorandum to the 
Applied Physics tLaboratory/Johns Hopkins University, 


Error Analysis When Brownian Motion is Used as a Datum 
Aging Model, by Barry BelKin, Jan. 13, 1978. 


wr a 
; Z 








We. 


ee 


is. 


14. 


Daniel H. Wagner, Associates, [nAtsrim Memorandum to 
Applied Physics Laboratory/Jonhns Hopkins University, The 
Ornstein-Uhlenbeck Displacement Process as a Model of 
Target Motion, by Barry BelKin, Feb. 1, 1976. 

Ibid. D. Si 

Ibid, p.22. 


Freedman, David A., Brownian Motion and Diffusion, p. I, 
Holden-Day, 19771. 


114 





INITIAL DISTRIBUTION LIST 


Defense Technical Information Center 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 6142 
Naval Postoraduate School 
Monterey, California 93946 


Department Chairman, Code 35 
Department of Operations Research 
Naval Postgraduate School 
Monterey, California 93948 


Dr. R. Neagle Forrest, Code 71 
Naval Postgraduate School 
Monterey, California 9%3946 


Dr. Donald P. Gaver, Code S5qQv 
Naval Postgraduate School 
Monterey, California 939498 


Dr. James N. Eagle, Code SSer 
Naval Postgraduate School 
Monterey, California 939406 


Dr. Kenneth V. Saunders 

Code 28231 

Strategic Systems Project Office 
Washington, D.C. 28376 


Mr. John C. Sommerer 

Submarine Technology Systems Group 
Applied Physics Laboratory 

Johns Hopkins University 

Johns Hopkins Rd. 

Laurel, Maryland 26767 


LT W. Justin Comstock 

c/o Dr. Sidney Q. Graves 
119 Dana Rd. 

Natchez, Mississippi 391280 


Copies 


bo 














Thesis 





omen 


200754 


C667 Comstock 


Gel 


Iwo models of time 
constrained target 
travel between two end- 
points constructed by 
the application of 
Brownian motion and a 
random tour. 


get tr 


cui 


ii 


Sue 
D 





i 


